Quantized conductance at the Majorana phase transition in a disordered 

superconducting wire 
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Superconducting wires without time-reversal and spin-rotation symmetries can be driven into a 
topological phase that supports Majorana bound states. Direct detection of these zero-energy states 
is complicated by the proliferation of low-lying excitations in a disordered multi-mode wire. We show 
that the phase transition itself is signaled by a quantized thermal conductance and electrical shot 
noise power, irrespective of the degree of disorder. In a ring geometry, the phase transition is signaled 
by a period doubling of the magnetoconductance oscillations. These signatures directly follow from 
the identification of the sign of the determinant of the reflection matrix as a topological quantum 
number. 
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It has been predicted theoretically [l[ that the s-wave 
proximity effect of a superconducting substrate can drive 
a spin-polarized and spin-orbit coupled semiconductor 
nanowire into a topological phase 0t3]) with a Majorana 
fermion trapped at each end of the wire. There exists 
now a variety of proposals [H-0) for topological quantum 
computing in nanowires that hope to benefit from the 
long coherence time expected for Majorana fermions. A 
superconducting proximity effect in InAs wires (which 
have the required strong spin-orbit coupling) has already 
been demonstrated in zero magnetic field [8] , and now the 
experimental challenge is to drive the system through the 
Majorana phase transition in a parallel field. 

Proposals to detect the topological phase have focused 
on the detection of the Majorana bound states at the end 
points of the wire, through their effect on the current- 
voltage characteristic [H, Gi| or the AC Josephson effect 
[ill [l2j . These signatures of the topological phase would 
stand out in a clean single-mode wire, but the multi- 
ple modes and potential fluctuations in a realistic system 
are expected to produce a chain of coupled Majorana's 
[l3l [l4 |. which would form a band of low- lying excita- 
tions that would be difficult to distinguish from ordinary 
fermionic bound states [l5j . 

Here we propose an altogether different detection strat- 
egy: Rather than trying to detect the Majorana bound 
states inside the topological phase, we propose to detect 
the phase transition itself. A topological phase transition 
is characterized by a change in the topological quantum 
number Q. The value of Q = ( — 1)'™ is determined by 
the parity of the number m of Majorana bound states 
at each end of the wire, with Q = — 1 in the topological 
phase jig . 

In accord with earlier work pj} , we relate the topolog- 
ical quantum number to the determinant of the matrix r 
of quasiparticle reflection amplitudes, which crosses zero 
at the phase transition. This immediately implies a unit 
transmission eigenvalue at the transition. Disorder may 
shift the position of the transition but it cannot affect 
the unit height of the transmission peak. We propose 
experiments to measure the transmission peak in both 



thermal and electrical transport properties, and support 
our analytical predictions by computer simulations. 

We consider a two-terminal transport geometry, con- 
sisting of a disordered superconducting wire of length L, 
connected by clean normal-metal leads to reservoirs in 
thermal equilibrium (temperature to). The leads sup- 
port 2N right-moving modes and 2N left-moving modes 
at the Fermi level, with mode amplitudes ip+ and ip—, re- 
spectively. The spin degree of freedom is included in the 
number N, while the factor of two counts the electron 
and hole degree of freedom. 

The 47V x 4N unitary scattering matrix S relates in- 
coming and outgoing mode amplitudes, 
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where the labels L and R distinguish modes in the left 
and right lead. The four blocks of S define the 2N x 2iV 



reflection matrices 



and transmission matrices t, t' 



Time-reversal symmetry and spin-rotation symmetry 
are broken in the superconductor, but electron-hole sym- 
metry remains. At the Fermi energy electron-hole sym- 
metry implies that if (u, v) is an electron-hole eigenstate, 
then also (v*,u*). Using this symmetry we can choose 
a basis such that all modes have purely real amplitudes. 
In this socalled Majorana basis S is a real orthogonal 
matrix, S t = <S* = S^ 1 . (The superscript t indicates the 
transpose of a matrix.) More specifically, since det<S = 1 
the scattering matrix is an element of the special orthog- 
onal group SO(4A). This is symmetry class D fl8l423j . 

The scattering matrix in class D has the polar decom- 
position 
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in terms of four orthogonal matrices O p G 
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SO (2 A) and 



a diagonal real matrix A with diagonal elements A„ G 
(— oo, oo). The absolute value |A n | is called a Lyapunov 
exponent, related to the transmission eigenvalue T n G 
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[0, 1] by T n = 1/ cosh 2 An. We identify 



2iV 



Q = sign Q, Q = Det r = Det r' = TT tanh A„ . (3) 



This relation expresses the fact that reflection from a Ma- 
jorana bound state contributes a scattering phase shift of 
it, so a phase factor of —1. The sign of n„tanhA n thus 
equals the parity of the number m of Majorana bound 
states at one end of the wire |24j . (It makes no difference 
which end, and indeed r and r' give the same Q.) 

To put this expression for Q into context, we first note 
that it may be written equivalently as Q = DetOiC?3 
if we restrict the A„'s to non- negative values and allow 
Det O p to equal either +1 or — 1. The sign of Q then 
corresponds to the topological classification of a class-D 
network model derived by Merz and Chalker [l 7 . We 
also note that Q can be written equivalently in terms of 
the Pfaffian of In M. Ai^ (with M. the transfer matrix in a 
suitable basis) [24( . A Pfaffian relation for the topological 
quantum number Q c iean in class D has been derived by 
Kitaev Q for a clean, translationally invariant system. 
We will verify later on that Q and Qcican agree for a 
clean system. 

An immediate consequence of Eq. ^ is that at the 
topological phase transition one of the A n 's vanishes 
[l7l [20L l2lj . so the corresponding transmission eigen- 
value T n — 1 at the transition point. The sign change 
of Q ensures that T n fully reaches its maximal value of 
unity, it cannot stop short of it without introducing a 
discontinuity in Q. Generically there will be only a sin- 
gle unit transmission eigenvalue at the transition, the 
others being exponentially suppressed by the supercon- 
ducting gap. The thermal conductance Gth = Go J2n Tn 
of the wire will then show a peak of quantized height 
Go = Tr 2 kgTo/6h at the transition. 

Our claim of a quantized conductance at the transi- 
tion point is consistent with earlier work [l9l - l22T ] on class 
D ensembles. There a broad distribution of the conduc- 
tance was found in the large- L limit, but the key differ- 
ence is that we are considering a single disordered sample 
of finite length, and the value of the control parameter 
at which the conductance is quantized is sample specific. 
We will now demonstrate how the peak of quantized con- 
ductance arises, first for a simple analytically solvable 
model, then for a more complete microscopic Hamilto- 
nian that we solve numerically. 

The analytically solvable model is the effective low- 
energy Hamiltonian of a class-D superconductor with a 
random gap, which for a single mode in the Majorana 
basis has the form 



H 



-ihvp<j z d / dx + A(x)a y 



(4) 



We have assumed, for simplicity, that right-movers and 
left-movers have the same velocity vf, but otherwise 
this is the generic form to linear order in momentum, 
constrained by the electron-hole symmetry requirement 
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FIG. 1: Thermal conductance and determinant of reflection 
matrix of a disordered multimode superconducting wire as 
a function of Fermi energy. The curves are calculated nu- 
merically from the Hamiltonian ©-(O on a square lattice 
(lattice constant a = i so /20), for parameter values W = l so , 
L = 10/so, A = 10-Eso, g c sf^BB — 21 E so , and three different 
disorder strengths Uo- The arrows indicate the expected posi- 
tion of the topological phase transition in an infinite clean wire 
(Uo = 0, L — ► 00), calculated from Eq. (|10|) . Disorder reduces 
the topologically nontrivial interval (where Detr < 0), and 
may even remove it completely, but the conductance quan- 
tization remains unaffected as long as the phase transition 
persists. 



H = —H*. An cigenstate $ of H at energy zero satisfies 



*(x) = exp 



Hvf 



A(x')dx' *(0). (5) 



By substituting *(0) = (I, 
the reflection amplitude 



•), *(X) = (i,0) we obtain 



r = tanh(LA//w F ), A = U 



A(x)dx. (6) 



In this simple model, a change of sign of the spatially 
averaged g ap A is the signature of a topological phase 
transition [251 ] . 

If A is varied by some external control parameter, the 
thermal conductance Gth = Go cos1i -2 (LA/?iuf) has a 
peak at the transition point A = 0, of height Go and 
width hvp/L (Thouless energy). The 1/ cosh line shape 
is the same as for a thermally broadened tunneling reso- 
nance, but the quantized peak height (irrespective of any 
asymmetry in the coupling to the left and right lead) is 
highly distinctive. 

For a more realistic microscopic description of the 
quantized conductance peak, we have performed a nu- 
merical simulation of the model |jj of a semiconduc- 
tor nanowire on a superconducting substrate. The 
Bogoliubov-De Gennes Hamiltonian 



H = 



A* E F 



A 



(7) 
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couples electron and hole excitations near the Fermi en- 
ergy Ep through an s-wave superconducting order pa- 
rameter A. Electron- hole symmetry is expressed by 

OyTyWOyTy = -% (8) 

where the Pauli matrices a y and r y act, respectively, on 
the spin and the electron-hole degree of freedom. The 
excitations are confined to a wire of width W and length 
L in the x—y plane of the semiconductor surface inversion 
layer, where their dynamics is governed by the Rashba 
Hamiltonian 

2 

H R = h U(r) + -^-{(JxPy - (TyPx) + \g c ttHBB(T x . 

2m e ff n 

The spin is coupled to the momentum p = —ihd/dr 
by the Rashba effect, and polarized through the Zee- 
man effect by a magnetic field B parallel to the wire 
(in the x-direction). Characteristic length and energy 
scales are l so — ti 2 /m c ga so and E so = m c sa 2 /h 2 . Typ- 
ical values in InAs are l so = 100 nm, E so = 0.1 meV, 
g e sHB = 2 meV/T. 

We have solved the scattering problem numerically [26[ 
by discretizing the Hamiltonian ([7|) on a square lattice 
(lattice constant a), with a short-range electrostatic dis- 
order potential U{x, y) that varies randomly from site 
to site, distributed uniformly in the interval (— Uo,Uq). 
(Eq uivalent results are obtained for long-range disorder 
(24|.) The disordered superconducting wire (S) is con- 
nected at the two ends to clean metal leads (Ni,N2), 
obtained by setting U = 0, A = for x < 0, x > L. Re- 
sults for the thermal conductance and topological quan- 
tum number are shown in Fig. U as a function of the 
Fermi energy (corresponding to a variation in gate volt- 
age). For the parameters listed in the caption the num- 
ber N of modes in the normal leads increases from 1 to 
2 at E F /E so w 10 and from 2 to 3 at E F /E so » 15. 
We emphasize that Fig. Q] shows raw data, without any 
averaging over disorder. 

For a clean system (Uq = 0, black curves) the results 
are entirely as expected: A topologically nontrivial phase 
(with Det r < 0) may appear for odd N while there is no 
topological phase for N even (27l - [29| . The topological 
quantum number of an infinitely long clean wire (when 
the component p x of momentum along the wire is a good 
quantum number) can be calculated from the Hamilto- 
nian rl(p x ) using Kitaev's Pfaffian formula 0, Hgj], 

Q clean = sign (Pf [a y T y H(0)]Pi [a y T y H(w/a)]) . (10) 

(The multiplication by o v T y ensures that the Pfaffian is 
calculated of an antisymmetric matrix.) The arrows in 
Fig. [T] indicate where Qclean changes sign, in good agree- 
ment with the sign change of Q calculated from Eq. 
(The agreement is not exact because L is finite.) 

Upon adding disorder Qclean can no longer be used 
(because p x is no longer conserved), and we rely on a 
sign change of Q to locate the topological phase transi- 
tion. Fig. [T] shows that disorder moves the peaks closer 




FIG. 2: Fourier amplitude with flux periodicity h/e of 
the magnetoconductance oscillations, calculated numerically 
from the Hamiltonian (O-© for a single disorder strength 
Uo — 50 -E S o and seven different temperatures To- The in- 
set shows the Aharonov-Bohm ring geometry. The parame- 
ters of the superconducting segment of the ring (S) are the 
same as in Fig. [1] with iV = 1 in this range of Fermi ener- 
gies. The normal part of the ring has N — 8 propagating 
modes to avoid localization by the disorder (which has the 
same strength throughout the ring). 

together, until they merge and the topological phase dis- 
appears for sufficiently strong disorder. We have also ob- 
served the inverse process, a disorder-induced splitting 
of a peak and the appearance of a topological phase, in 
a different parameter regime than shown in Fig. [TJ Our 
key point is that, as long as the phase transition persists, 
disorder has no effect on the height of the conductance 
peak, which remains precisely quantized — without any 
finite-size effects. 

Since electrical conduction is somewhat easier to mea- 
sure than thermal conduction, we now discuss two al- 
ternative signatures of the topological phase transition 
which are purely electrical. An electrical current Zj is 
injected into the superconducting wire from the normal 
metal contact Ni, which is at a voltage V\ relative to 
the grounded superconductor. An electrical current I2 
is transmitted as quasiparticles into the grounded con- 
tact N2, the difference I\ — I2 being drained to ground 
as Cooper pairs via the superconductor. The nonlocal 
conductance G — I2/V1 is determined by the time aver- 
aged current I2, while the correlator of the time depen- 
dent fluctuations SI2 determines the shot noise power 
P = f^dtiSh^Shit)) (in the regime k B r < eV 1 
where thermal noise can be neglected). 

These two electrical transport properties are given in 
terms of the N x N transmission matrices t ee and the 
(from electron to electron and from electron to hole) by 
the expressions [3(j 

G=(e 2 /h)TrT-, P = (e 3 Vi/h) Tr (7+ - T 2 ), (11) 
r±=tl e t ee ±t{ e t he . (12) 

Electron- hole symmetry relates t ee = t* hh and the = t* h . 
This directly implies that Tr7+ = ±Trttt = \ Y, n T n- 
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If in addition we assume that at most one of the T n 's is 
nonzero we find that 71 vanishes [24] . We conclude that 
G remains zero across the topological phase transition, 
while P/Vi peaks at the quantized value e 3 /2h. This is 
the second signature of the phase transition |3l| . 

The third signature is in the electrical conductance. 
Since G = for a single open transmission channel, we 
add (topologically trivial) open channels by means of a 
parallel normal metal conductor in a ring geometry. A 
magnetic flux $ through the ring produces Aharonov- 
Bohm oscillations with a periodicity A$ = h/e*. The 
effective charge e* = e if electrons or holes can be trans- 
mitted individually through the superconducting arm of 
the ring, while e* = 2e if only Cooper pairs can be trans- 
mitted [HI, EH . We thus expect a period doubling from 
h/2e to h/e of the magnetoconductance oscillations at 
the phase transition, which is indeed observed in the 
computer simulations (Fig. [5J . To show the relative ro- 
bustness of the effect to thermal averaging, we repeated 
the calculation at several different temperatures tq. For 



E so ~ 0.1 meV the characteristic peak at the phase tran- 
sition remains visible for temperatures in the readily ac- 
cessible range of 100-500 mK. 

In conclusion, our analytical considerations and nu- 
merical simulations of a model Hamiltonian [l| of a dis- 
ordered InAs wire on a superconducting substrate show 
three signatures of the transition into the topological 
phase (Figs. [T] and [2]): A quantized thermal conductance 
and electrical shot noise [3l|, and a period doubling of 
the magnetoconductance oscillations. These unique sig- 
natures of the Majorana phase transition provide alter- 
natives to the detection of Majorana bound states [§- 
EH, EH, which are fundamentally insensitive to the ob- 
scuring effects of disorder in a multimode wire. 
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Appendix A: Derivation of the scattering formula 
for the topological quantum number 

1. PfafHan form of the topological quantum 
number 



The topological quantum number Q of a disordered 
wire is given in Eq. ([3]) as the sign of the determinant of 
the reflection matrix. That is the form which is most con- 
venient for computations. In order to derive this relation- 
ship and also to compare it with results in the literature 
for translationally invariant systems it is convenient 
to rewrite it in terms of the transfer matrix M. It then 
takes the form of a Pfaffian, rather than a determinantal, 
relation. 

The 4A x AN transfer matrix M relates the mode am- 
plitudes to the right (R) and to the left (L) of the disor- 
dered wire, 



E 



= M 



V>-,L 



M = 



(Al) 



The condition of particle current conservation is 
a z M^a z = M^ 1 , where the Pauli matrix a z acts on the 
block structure indicated in Eq. (|A1[) . In the Majorana 
basis of real mode amplitudes M is a real matrix, hence 

a z M t cr z = MT 1 . (A2) 

The transfer matrix has the polar decomposition 



M = 



o 2 



cosh A — sinh A 
- sinh A cosh A 



3 

o o\ 



O2 \ , . > fo 3 



(A3) 



where the matrices O p G SO (2 A) and A = 
diag (Ai , A2 • • • , \2n) are the same as in the polar decom- 
position ^) for the scattering matrix. One readily checks 
that Eq. (|A2I) is satisfied. 

The polar decomposition (| A3|) allows us to compute 



n = 



o 2 



This is an antisymmetric matrix, so it has a Pfafhan 

2JV 



Pf (a z \nMM^) = JJ 2A.,, 



(A4) 
a, 

(A5) 



We have used the identity 

Pf BAB* = Det B Pf A, 



(A6) 



with Det CI = 1. 

We conclude that the topological quantum number ^ 
can equivalently be written as 



Q = sign [Pf (a z In MM*)] 



E 




(A7) 



' ordinary fermion 
Majorana fermion 



FIG. 3: Procedure to count weakly coupled Majorana bound 
states in a disordered superconducting wire. Majorana 
fermions at the two ends of the wire (light blue) are weakly 
coupled, so their energy is not exactly zero and we need a 
way to distinguish them from an ordinary fermionic excitation 
(dark blue) . To that end we close the wire into a ring through 
a topologically trivial superconductor and ask whether de- 
structive interference of the tunnel splitting in the two arms 
of the ring can produce a pair of two-fold degenerate zero- 
energy states. 



2. How to count Majorana bound states 

To determine the topological quantum number of the 
disordered superconducting wire we seek the number 
of Majorana bound states. Particle-hole symmetry en- 
sures that any bound state at zero energy is a Majorana 
fermion (since the creation and annihilation operators are 
related by 'y'(E) = j(—E) and therefore are identical at 
E = 0). However, we cannot directly search for zero- 
energy eigenstates: Even if the Majorana fermions are 
maximally separated by the entire length L of the wire 
they will still have a nonzero tunnel coupling which splits 
their energies apart, away from zero. 

The issue here is how to distinguish strongly coupled 
from weakly coupled Majorana fermions. Any ordinary 
fermionic excitation, with distinct creation and annihi- 
lation operators a* 7^ a, can be described by a pair of 
strongly coupled Majorana fermion operators 71 = a + a\ 
72 = i(o. — a)). In contrast, the Majorana bound states 
at opposite ends of the wire are weakly coupled Majorana 
fermions. 

Our geometry of a disordered wire connected at the 
ends to metal contacts allows for a natural distinction of 
weak versus strong coupling: We call a pair of Majorana 
bound states "strongly coupled" if they are more strongly 
coupled to each other than to one of the ends of the wire. 
Conversely, weakly coupled Majorana bound states are 
more strongly coupled to one end of the wire than to 
any other Majorana. The topological quantum number 
counts only weakly coupled Majorana's. 

This distinction between weak and strong coupling can 
be made operational by means of the thought experiment 
illustrated in Fig. [3j We close the wire into a ring by 
connecting the two ends through a superconductor which 
is in a topologically trivial phase (with a uniform positive 
gap An). Destructive interference in the two arms of the 
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ring can eliminate the tunnel splitting between a pair of 
Majorana bound states and produce two-fold degenerate 
zero-energy eigenstates, if the coupling between the two 
Majorana's through each arm of the ring is of comparable 
strength. 

So we vary Ao (allowing also for mode mixing at the 
junction between the two arms of the ring) and find that 
a number m of two-fold degenerate states appear at zero 
energy. This means that the disordered wire contains 
m pairs of Majorana's which are more strongly coupled 
to the ends of the wire than to each other (otherwise 
the couplings through the two arms of the ring could not 
have been equalized by varying Ao). The number m thus 
counts the number of weakly coupled Majorana bound 
states, which gives the topological quantum number Q = 
(-l) m 



3. Topological quantum number of a disordered 
wire 

Now that we have an operational definition of the topo- 
logical quantum number of a finite system, our next step 
is to relate this to the scattering parameters A„ in Eq. 
(U). For this purpose it is easiest to work with the trans- 
fer matrix, rather than the scattering matrix. An eigen- 
state of the ring must be single- valued as we go around 
the ring, so in terms of the transfer matrices A4 and 
M.q of the two arms of the ring we have the condition 
AdoM^ — 'F This leads to the determinantal condition 



Det(l-M M) =0. 



(A8) 



We choose to work in a basis where the orthogonal 
matrices O p in Eq. (|A3I) are equal to the unit matrix. 
Each of the n = 1,2,... 2N eigenchannels of the disor- 
dered wire can then be treated separately, with 2x2 
transfer matrices M n = exp(— cr x A n ) at zero energy. The 
topologically trivial arm of the ring (of length Lg and co- 
herence length £o = Hvf/Aq > 0, without any disorder) 
has transfer matrix Mo = cxp(— cr x Lo/^o). The condition 
for an eigenstate at zero energy reads 



Det (1 - e ~°* L o/io e -<?*^ = 0) 



(A9) 



which has a twofold degenerate solution if the ratio Lq / Co 
is tuned to the value — A„. This is the pair of weakly 
coupled Majorana bound states in the n-th eigenchannel 
that we are searching for. Because £o > 0, by definition 
in a topologically trivial phase, the pair exists only if 
K < 0. 

We conclude that the number of pairs m of weakly cou- 
pled Majorana bound states equals the number of nega- 
tive A„'s, hence 



Q = (-ir = sign (ll A «) 



(A10) 
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FIG. 4: Same as Fig. [U for an impurity potential with cor- 
relation length £ = 2 a (upper panel) and £ = 10 a (lower 
panel. 



Appendix B: Numerical simulations for long-range 
disorder 

Fig.CQinthe main text demonstrates that the quantized 
thermal conductance at the Majorana phase transition is 
insensitive to short-range disorder (correlation length £ 
of the order of the lattice constant a) . Here we show that 
long-range disorder similarly has no effect on the quanti- 
zation. (The stability of Majorana bound states against 
short-range and long-range disorder was investigated in 
Ref. [1|.) 

As before, we solve the scattering problem numerically 
by discretizing the Hamiltonian (J7J on a square lattice 
(with a total number of N tot lattice points in the disor- 
dered region). The disorder is modeled as a superposition 
of impurities with a Gaussian profile, 



Wimp 

U(r) = ^ ex P 

i=l 



(r 



2? 



(Bl) 



as announced in Eq. ([3]). 



where -/V; mp is the number of impurities. (We fixed the 
impurity concentration 7ij mp = iVj mp /iVtot at 5%.) The 
strength U% of an individual impurity is randomly dis- 
tributed in the interval (—Uq,Uo), and the impurity posi- 
tions Ti are chosen randomly from the iVtot lattice points. 

The results of the calculation are shown in Fig. [4j for 
different values of the correlation length £. In all cases, 
we observe as before that the thermal conductance re- 
mains quantized as long as the topological phase persists. 
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For sufficiently strong disorder, the merging of two peaks 
signals the disappearance of the topological phase and a 
breakdown of the conductance quantization. 

Appendix C: Electrical conductance and shot noise 
at the topological phase transition 

The expression (ITT|) for the nonlocal electrical conduc- 
tance and shot noise of the superconducting wire can be 
evaluated further if there is only a single open transmis- 
sion channel. The 2N x 2N transmission matrix 

t= (\ ee i eh \ (ci) 

\tfie thh J 



is then of rank 1, which means that the N x N subma- 
trices t ee , thh, the, teh have the dyadic form 

tee = |ur)(«l|, t h h = |«r)(«l|i 

the = |«r)(«l|, t e h = \un)(vh\- (C2) 

The matrix 7± then becomes 

T± = |ul}<ul|((ur|ur) ± (va|va)). (C3) 

Electron- hole symmetry requires \vn) = \u^), hence 
71 = 0, T+ = iTrttt ; anc i thus G = 0, P = 
(e 3 i/i/2/t)Trttt. 



